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Analytical solutions to the chaotic and ergodic motion of a certain class of one-dimensional 
dissipative and discrete dynamical systems are derived. This allows us to obtain exact expressions 
for physical properties like the time correlation function. We illustrate our solutions by means of 
a few examples for which conventional numerical trajectory calculations fail to predict the correct 
behaviour. 
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I. INTRODUCTION 



Chaotic behaviour is a characteristic feature of the overwhelming majority of deterministic nonlinear dynamical 
systems pi. Classical chaos manifests itself in the exponential sensitivity of the trajectories with respect to the initial 
^ ■ conditions and has nicely been understood as the very complex topological process of ever lasting stretching and 
Cn I folding of the motion. In view of the complexity of chaotic trajectories it appears hard to imagine that chaotic motion 
'— ' ' could be described by means of closed analytical formulae. To our knowledge there exists only one exception, i.e. 
I~C ' dynamical law, which allows the representation of its chaotic trajectories for arbitrarily long times in a closed analytical 
I form: the solutions of the logistic map in the ergodic chaotic limit are given by the celebrated Pincherle relation [E| . 

h — In contrast to our lack of knowledge concerning exact expressions for chaotic motion there exists a well-founded and 
^\ : justified interest in obtaining solutions to ergodic and chaotic behaviour. 

On the o.h.s. it is clear that the frequently used numerical solutions of chaotic dynamical systems do not yield one 

and the same trajectory for long time scales S. Nevertheless due to the shadowing argument statistical quantitites 

i-rt '' of chaotic ergodic systems with and without external noise can in many cases approximately be obtained through a 

J^ ' numerical investigation [EUs] . However there are exceptions like the case when a Lyapunov exponent fluctuates about 

zero and this is expected to be common in simulations of higher-dimensional systems B . It is therefore highly desirable 

fH to make exact properties accessible: if the exact trajectory could be derived this would be an excellent starting-point for 

^ , calculating the exact correlation function, invariant measure, Lyapunov exponent or other quantities [Q . In particular 

IL/ ' it would also offer the possibility of determining the exact long time behaviour of relevant physical quantities. Even 

• I— I , more important may be the fact that analytical representations of chaotic motion are of principle interest and can 

J^ ' certainly enhance our understanding of the complexity of chaotic dynamics in general. Typical questions which could 

JH ] then be addressed are: How does the exponential sensitivity with respect to the initial conditions come about? Are 

. - - 1 there any characteristic (scaling) properties and/or self-similar structures of the analytical expressions which are 

responsible for the complex stretching and folding process of chaotic dynamics? 

The current investigation represents a first step in the above direction and provides closed analytical formulae for 
two classes of one dimensional unimodal, dissipative, ergodic and chaotic maps of the interval each containing an 
infinite number of members. In particular we will derive exact solutions for two classes of maps: the conjugates 
to both the symmetric as well as asymmetric tent maps. Using these solutions we will derive exact expressions for 
important physical quantitites like, for example, the correlation function. Some specific examples are discussed in 
detail thereby demonstrating the deviation of exact and numerically determined properties. 

II. SOLUTIONS TO THE CHAOTIC DYNAMICS OF THE GENERAL CONJUGATES OF THE 

SYMMETRIC TENT MAP 

A. Trajectories 

Let us begin our investigation by considering the symmetric tent map (STM) 0. The n-th iterate of the STM 
t'^^\x) possesses 2" monotonicity intervals with alternating constant slopes 2" and —2". The zeros and maximas are 



at the positions ^ with k = 1, ..., 2". This makes it possible to represent t'^^\x) in the form 

The above formula can be recast into the following very simple expression 

t(") (x) = 1 - I (2" a; mod 2) - 1 1 (2) 

where mod is the modulo operation. We consider now the family of maps g{x) = uot o u^^{x) which are obtained 
from the STM by conjugation with an invertible and differentiablc function u{x) which maps the unit interval onto 
itself and obeys w(0) = 0, m(1) = 1 ||]. The rt-th iterates of the maps {g{x)} are then given by 

g(")(x) — (uoto u^^) o [uoto u^^) o...o[uoto u^'^){x) = uo i^") o u^^{x) (3) 

and therefore 5^"-' (a;) takes on the following appearance 

g^"'>{x)=u{l-\{{2"u~^{x)) mod2)-l|) (4) 

which represents a closed form analytical solution to the dynamics of the maps conjugated to the STM. The general 
conjugation yields a variety of symmetric as well as nonsymmetric maps. Imposing the restrictive condition u(x) — 
1 — u{l — x) for the conjugating functions u{x) we obtain the special class of so-called doubly symmetric maps for 
which both the invariant density as well as the map is symmetric. In particular if we use the specific conjugation 
ui{x) = sin^ (^) and correspondingly u^^{x) = i arccos(l — 2x) we arrive at the logistic map l{x) = 4a;(l — x) 
and after a little algebra at the Pincherle relation l^"''{x) = ^ {1 — cos [2" arccos(l — 2a;)]} for the n-th iterate of the 
logistic map. 

Eq.(4) nicely demonstrates the exponential sensitivity with respect to the initial conditions as well as the infinite 
process of stretching and folding. The stretching process takes place through the multiplication with an exponential 
factor (2") and the folding through the mod- function which cuts with increasing iterations an increasing number of 
digits making a more detailed specification of the initial conditions necessary in order to describe the actual motion. 

B. Exact Properties and Examples 



Since the Lyapunov exponent is invariant with respect to conjugation [BUgf| all members of the above-discussed class 
of maps possess the same Lyapunov exponent A — In 2. In particular, since /x{a;) = a; is the invariant measure of the 
STM, the measure of the conjugated maps is fi{x) = u^^{x) and, therefore, varies widely with changing conjugating 

function. Next let us derive analytical expressions for another important physical quantity namely the correlation 

1 1 

fiuiction which is defined by C{n) — C{n) — xP' with C{n) — J g^"'{x)x dfi{x) and the mean value x = J xdfi{x) Let 



us denote the positions of the zeros and maxima of g^^'{x) by {x2k} and {a;2fc-i}, respectively. Using the fact that 

g^"'>{x) is conjugate to some 'original' map h^"\x) as well as the relation u{yi) — Xi the correlation function can be 

decomposed as follows 



J/2fc + l V2{k+1) 

C{n) ^ y^^ Ck{n) with Ck{n)^ f u{x)u{h^"\x)) dx + f u{x)u{h^"\x)) dx n>0 (5) 

y2k V2k+i 
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In order to provide exact correlation functions for some specific classes of maps we now specialize to the case that 
h{x) is given by the STM. Using eq.(l) we can derive the following simple structure for the terms of the sum of the 
correlation function 

2fc-l k 

2" 2^-1 

Ck{n)^ I u{x)u{2''x-2{k-l))dx+ [ u{x)u{~2"x + 2k)dx (6) 



In general the terms Ck{n) can be evaluated analytically as we shall show in the following by means of a few 
examples. Let us first choose the conjugation u{x) — a;i+^ which results in the invariant measure /i(a;) — x^^^ and 



the corresponding normalized power law density p{x) — (/? + l)x^ , (3 > —1. Due to their scaling properties power 
law densities are of particular interest for physical systems with critical or self-similar behaviour. Using eqs.(5,6) we 
arrive after some algebra at the following closed form analytical expressions for the corresponding correlation function 



(3 + l\ f l\^ , ,^-:(^+^) , ,(i_„)4±? V- fl.. .N^oA /3 + 2 
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where B{a,b) and B{a,b,x) denote the complete and incomplete unnormalized Beta function, respectively, and 
2Fi{a,b,c,x) is the hypergeometric function. For the particular case P — the above expression reduces to the 
(^-correlation C{n) oc (5„o as expected. Figure la shows the absolute value of the above correlation function for differ- 
ent values of the power /?. The short time as well as long time behaviour changes somewhat with changing parameter 
/?. The asymptotic behaviour (n — > oo) of C(n) is an exponential decay. The decay constant r can be determined 
using the Euler-MacLaurin sum formula |10|| for the asymptotic expansion of eq.(7). It turns out that for /3 > 
T = 4^1 In 2 while for — 1 < /? < we have r — 2 In 2 independent of /?. Having obtained the exact correlation 
function for the above class of maps with power law density we are in the position to compare these results with those 
of numerical trajectory simulations. Figure lb shows the comparison of the exact correlation function for (3 = 5.0 (full 
circles) with the results of numerical trajectory calculations using 10^ (squares) 10'' (diamonds) and 10^ (triangles) 
points. Obviously there is a strong inherent deviation and in particular it can be observed that an enhancement of the 
number of points of the numerically calculated trajectory does not yield an improvement in the statistical accuracy of 
the corresponding correlation functions. This is due to the fact that a naive numerical simulation looses for the above 
dynamical systems very rapidly the original trajectory and gets due to numerical inaccuracies trapped on a certain 
'orbit'. 

Finally let us provide one more example out a large number accessible by the above given analytical formulae. We 
choose as a conjugating function u{x) ~ sin(^x) which results in the measure /i(a;) — - arcsin(x) and the invariant 
density p{x) — — ,}' 2 ■ ^ straightforward but tedious calculation yields the terms Ck{n) which can be summed up 
to the final beautiful result of the correlation fmiction 

C{n) = ,f^' , cot (^) - (-) (8) 

which shows with respect to its asymptotic behaviour (n -^ 00) an exponential decay with a decay constant 2 In 2. 

III. SOLUTIONS TO THE CHAOTIC DYNAMICS OF THE CONJUGATES OF THE ASYMMETRIC 

TENT MAP 

Solving the problem of the chaotic and ergodic dynamics generated by the general class of maps conjugated to the 
asymmetric tent maps (ATM) S is a much more intricate task than the corresponding solution of the dynamics of 
the maps conjugate to the STM. In order to develop the necessary concepts and techniques we therefore proceed in 
several steps. 

A. The Concept of Mode Superposition 

The knowledge of the position of the minima/maxima, i.e. the length of the monotonicity intervals, of the n-th 
iterate of the ATM is central to the derivation of an analytical expression for the n-th iterate of the ATM. In the 
present subsection we therefore provide a method which allows to determine the length Ln(i) of the i-th monotonicity 
interval (we count the intervals starting from the origin) by a decomposition technique of the corresponding frequency 
distribution r„(«) (see below) into independent modes. 

The ATM is unimodal with constant slopes 1/p and ~l/q {q — (1 — p)) on the two monotonicity intervals, re- 
spectively. First of all we remark that each iteration process n ^ n + I divides a monotonicity interval of the n-th 
iterate into two monotonicity intervals of the {n + 1) — th iterate which possesses positive and negative slopes in these 
intervals, respectively. The ratios of their lengths are given hy p : q or q : p depending on whether the n-th iterate of 
the ATM on the original monotonicity interval possesses a positive or negative slope, respectively. In order to specify 



the length of the i-th monotonicity interval we do not need such detailed information like the symbolic sequence |^ 
but only the frequency of the occurence of the factors p and q during the iteration process. The number of different 
lengths for the monotonicity intervals is therefore much smaller than the number of monotonicity intervals themselves. 
The length of the i-th monotonicity interval of the n-th iterate of the ATM is given by 

in(i)=p"-''"(')g''"(^^ (9) 

where r„ (i) is the frequency of the occurence of the factor q during the branching process into the final i-th monotonic- 
ity interval of the n-th iterate. Using the above-described properties of the branching of the monotonicity intervals it 
can be shown that r„(i) obeys the following recursion formula 

ri(l) = 

..^ /r„(z) ; l<z<2" 

'^"+1^*^" \r„(2"+i-(?:-l)) + l ; 2" + l<z<2"+i ^^^' 

r„+i(i) can therefore be obtained from r„(i) through a reflection cr„ of the function r„(i) (i — 1, ..., 2") with respect 
to the vertical axis located at 2" -I- ^ and subsequent addition t„ of 1 to the attached part of r„+i(i) which yields in 
total r„+i(i) = (r„ o cr„)r„(i). This process is illustrated in Figure 2 for the three functions ri(i), r2(i), ?'3(i). 

In the following we derive a decomposition of the functions ?"„(«) into different modes, i.e. r„(i) can then be 
described by a superposition of these modes. Let us introduce n — 1 modes Mn,m{i)i'm — 2. . .n and the additional 
mode Nn- Each mode Mn.m{i) is defined on the finite support i = 1, ..., 2" and oscillates as a function of i with the 
period Tm = 2'". In addition it possesses a phase (j>{Mn,m) = 2™^^ which describes the shifting of the oscillations 
on the i-axis. The mode Nn is characterized by the period 2" and the phase 0(-/V„) = 2"~^. The equally weighted 
superposition of these modes yields the quantity r„(i), i.e. we have 

n 

rn{^) = Yl ^'^".™(*) + ^"(0 (11) 

Figure 3 illustrates the three modes M3_2(i), ^3^3(1) and N3{i) which are necessary in order to build up r3(i). In 
particular we now briefly verify that the above decomposition obeys the recursion formula given in eq.(lO). The action 
of the reflection (T„ onto the modes Mn,m{i) is given by 

M„,„ ^^ Mn+i,m Nn -^ M„+i,„+i (12) 

The application of the subsequent addition t„ to the reflected modes corresponds to the inclusion of the mode Nn+i- 
In total we therefore arrive at the desired relation 

n n+1 

r„(i) = Yl Mn^t) + Nn{i) "^^" r„+i(z) = ^ M^+iM^) + Nn+i{i) (13) 

m—2 m—2 

In the following subsection we provide construction principles and closed analytical formulae for the modes Mn,m{i) 
and Nn- 

B. Construction of the Modes 

Each mode function is defined on 2" natural numbers and their continuation and representation on the real axis 
is therefore not unique. In the present subsection we provide two different representations for the mode functions 
Mn,m{i) and Nn{i)- The first one is characterized by the application of the step function to trigonometric functions 
and the second one uses polynomials and their periodic continuation in order to describe the modes. 

Looking at Figure 3 suggests a representation of the modes by the action of the step function {Q{x) — 
{1 : cc > ; : a; < 0}) onto an oscillating function. Choosing the sin function for the oscillating part we can 
adapt the periodicity and phase in order to obtain the modes M„ m(*) and Nn{i), respectively. They take on the 
following appearance 



^{^) = e 



2i - (2™-i + 1) 
vr 



Nn{i) = e 



2i - 2" + 1 

sm I TT ^ - 

2" 



(14) 



For the function r„(i) wc therefore arrive at the expression 



.W 



m=2 



e 



sm TT 



e 



_2i- (2" + l) 
2" 



(15) 



Introducing the above formula for r„(i) in eq.(9) we have an analytical expression for the length of the i-th monotonicity 
interval. The only feature which could be seen as a disadvantage of the above representation is its discontinuity through 
the introduction of the step function. It is therefore desirable to gain a second representation of the mode functions 
which should have the property of being smooth with respect to the continuation onto the entire real axis. We herefore 
use polynomials whose coefficients will be adapted correspondingly. The introduction of a sin function yields then a 
periodic continuation of the polynomials and allows to adapt to the functional form of the modes via a corresponding 
scaling operation and phase shift. 

The first step in the derivation of a second, smooth, representation of the modes Mn^m{i) and Nn(i) is the con- 
struction of a polynomial with suitable properties. To achieve this, our starting-point is a polynomial Pmix) which 
possesses at 2™"^ points xj' the value +1 and at other 2™^^ points x^ the value —1. 



Pm\^l ) — Pm\^2 ) — • • • — Pm\^2^-^ ' — Pm\^i ) — Pm\^2 ) — ' • " — Pm\^2^-^ ) — 



(16) 



For the time being the domain of definition of the above introduced polynomial includes 2™ integer values and 
possesses 2™^^ zeros. Such a polynomial can be written as a sum of 2"* terms, each of the terms thereby acquires the 
value -1-1 or —1 at exactly one position given by xf or x^ , respectively whereas at all other positions x'^ ,xj , i ^ j 
it vanishes. These properties are guaranteed if each term consists of a normalized product of linear factors and we 
therefore arrive at the following expression for the polynomial Pm{x) 



Prn{x) = ^ 



(x-x^) 



n 



{x-x+)ix-x^ ) 



{x 



n 



{x-x+){x-x. ) 



t i^t-^^) 7=1 i^i -^j)i^i -^j) ~t i^i -^i) 7"=i i^2 -^J)iXt -^j) 



(17) 



The points of support xf , x^ ought to be equidistant in the unit interval [0, 1] and should be arranged symmetrically 
with respect to 0, i.e. we have 



2i~l 



(18) 



After a little algebra this leads to an essential simplification of the polynomials (17) 



Pmix) = 



Z' 



^n 






As a next step we perform a periodic continuation of the polynomial pm{x) by substituting 

X > sin(7rx) Xi > sin(7rxi) 



(19) 



(20) 



In order to describe the modes with the constructed polynomials we have to adapt the frequency as well as phase of 
the oscillating periodic functions defined by eqs.(19,20). In addition a shift is performed in order to make the points 
of support equal to integers on which the modes are defined. Finally we arrive at the following expressions for the 
individual modes 



Mn,m{x) 



E 



.2a;-(2'""Hl) 



Sin {n^ 



n 

j=i 



2x-(2'""Hl) 



— sin^ (tt 



sm^f^^i^ 






sm 




(21) 



Nn{x) = 




sm n 



. 2a-(2'' + l) \ 2"- 



sin (tt^I^ 



2" 
2" 



n 

3 = 1 



1 ■ 2 ( 

sm I 71" 



2a:-(2"+l) 
2" 



sm 



2j-l ^ 



{-' 



{- 
(-^) 



(22) 



and for the frequency distribution r„ (i) of the i-th nionotonicity interval 

f 2-"- 3in (n^-%Z'^'A --^ Bin^ U-^^^K^^ sin^ (-^) 



^w-E^^+E \,,(.^. ^ n 



3 = 1 



sin^ (.21^) - sin^ {n^^) 



+ - 



y ^in (^^^^^i;^) ^- sin^ (^^^^^g;^) - sin^ (^^) ^ 
^'+^ sin(.^) fi sin^ (.21^) - sin^ (.2^) ^ 



This concludes our construction of the modes. 



C. Analytical Representation of the Chaotic Dynamics 



(23) 



The knowledge of the frequency distribution r„(i) is a key ingredient for the derivation of analytical solutions to 
the chaotic motion of the conjugates of the ATM. Arbitrarily high iterates are now accessible in closed analytical form 
via the following procedure. 

First of all we give some important quantities related to the n-th iterated map: 



The length of the i-th nionotonicity intervals is given by 



(24) 



• The slope in the i-th nionotonicity interval is 



S'„(i) =p'-"W-"g-'-"«(-l)(^-i) 



(25) 



• The left and right boundary G^(j) and G^{i) of the i-th nionotonicity interval can be obtained by summation 
of the lengths of the corresponding intervals: 



i-l 



G^(z) = ^L„(j)=5]p"-'-"(^')g'-"(-'") 



i=i 



i=i 






(26) 



• Definition of the index-function: Let us define a function /„ : [0, 1] — > {1, ..., 2"} whose value In{x) tells us the 
number of the nionotonicity interval in which x is contained. The index function can be described with the help 
of the step function Q{x) in the following way 



2" 



/„(x) = ^e[x-G^(z)]=5]e 



i=l 



i-l 



(27) 



We are now in the position of formulating the n-th iterate tp' (x) of the ATM in a closed analytical expression 

t(")(x) ^{x- G^iUx)]} S4Ux)] + 1 (l + (-1)^"(-)) (28) 

Inserting eqs. (26,27) into the above expression yields 

/„(x)-i 



t^i^) 



y^ p"-''"0')(,''"0') 



j=i 



^rr.ll^ix)]-n^-r„[l^ix)]^_^y„ix)-l _^_ }_ f i _^_ (_l)^-.(^)') 



(29) 



»t 



Similar to the case of the STM we obtain the n-th iterate gl (x) of the maps conjugated to the ATM by the following 
decomposition 



5t (a;) = u o tp o u o uotp . . . otp o u {x) — u o v^' o u {x) 



(30) 



which gives together with eq.(29) our final result 



at^^) = "' 



7„(«-i(x))-l 



s., 



(/„(u-i(x))) + i(l + (-l/"(-^(^))) 



(31) 



Eq.(31) can be used to calculate arbitrarily high iterates of the conjugates of the ATM. As an example we show in 
figure 4 the 8-th iterate of the conjugate u{x) = siv? {^) of the ATM for p = 0.9 which gives an idea of the scaling 
structures contained in the above analytical formula. 

The invariant measure of the conjugates of the ATM is given by fi{x) = u~^{x) and the Lyapunov exponent is 
known to be A = —pln{p) — (l—p) ln(l —p) [Q. The above eq.(31) can be used to calculate the correlation functions 
of the conjugates of the ATM as well as any other physical property. 

IV. CONCLUSIONS 

We have presented analytical solutions to the chaotic and ergodic dynamics of two classes of noninvertible single 
humped maps: the conjugates of the symmetric as well as asymmetric tent map. With the exception of the iterates of 
the logistic map this is, to our knowledge, the first time that closed analytical representations of chaotic trajectories 
are derived. This enables us to calculate exact physical quantities, like, for example, the time correlation function, 
even in those cases for which numerical trajectory calculations fail to predict the correct long time behaviour. The 
complex stretching and folding process of the iteration of the dynamical system is clearly revealed within our analytical 
approach through the combination of conjugating functions, scale transformations and cut-off operations. 

The key to the construction of the analytical solutions is the complex superposition of modes connected with the 
monotonicity intervals of the n-th iterate of the dynamical system. We conjecture that this holds also for the general 
case of arbitrary single humped maps and might be a hint towards the construction of their solutions. The boundaries 
of the monotonicity intervals are then given by the preimages of the maximum of the dynamical law. 



FIGURE CAPTIONS 

Figure 1: (a) The logarithm of the absolute value of the correlation function for the conjugation u{x) = a; i+^ of the 
STM for different values of the parameter /3. The values f3 = —0.95, —0.5, +1.0, -1-5.0, -1-50.0 correspond to the curves 
with the full circles, circles, plus, full triangles and full boxes, respectively, (b) The exact correlation function for 
f3 = 5.0 (solid line) together with the correlation functions resulting from numerical trajectory calculations (broken 
lines) with 10"^ (squares), lO** (diamonds) and 10^ (triangles) points. 

Figure 2: The frequency distributions ri{i),r2{i) and r^{i) of the ATM. 

Figure 3: The modes Af3,2(«): ^^3.3(*) and N^ii) of the ATM. 

Figure 4: The 8-th iterate of the conjugate u{x) — sin^ (^) of the ATM for p — 0.9 calculated through eq.(31). 
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